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It has been a central topic of surface science how to 
control the adsorption and desorption in order to to 
bring out desirable features and functionalities by ad- 
sorbed molecules. Tuning the electronic features of solid 
surfaces has an important implication in that molec- 
ular chemisorptions and catalytic reactions are deter- 
mined by them 1 ' 2 . In particular for graphene, the two- 
dimensional honeycomb carbon lattice, in which the con- 
duction 7r*-band and the valence 7r-band contact to each 
other at the "Dirac point" giving a feature of zero-gap 
semiconductor 3 , the control of chemisorption is a critical 
issue since chemisorption directly leads to altering every 
electronic property of graphene. Other than the elec- 
tron/hole doping 4 owing to the charge transfer between 
graphene and the adsorbed molecules, widely known are 
the charged impurity effect on the electron transport 5 7 , 
lattice deformation 8 , and opening the band gap due 
to asymmetric adsorption 9 11 . Aside from the macro- 
scopic spatially-controlled adsorption that is achieved us- 
ing nano-device fabrication technique 12 , microscopic con- 
trol of adsorption structure is of great importance be- 
cause the aforementioned adsorption effects are altered 
by the local structure of adsorbate, e.g., whether the 
adsorbed molecules are arranged in a random or su- 
perlattice structure 10 ' 13 , or whether the molecules are 
adsorbed individually 4 ' 14 or collectively (in dimers 15 or 
clusters 16 ' 17 ). For the first step to realize such an 
advanced control of adsorption, the methods to uti- 
lize the interaction between the adsorbed molecules and 
graphene for it are to be explored. 

The principal impetus in the present study is to con- 
trol the charge transfer between graphene and the ad- 
sorbed molecules by tuning the Fermi level of graphene, 
which is readily accomplished in the field effect transis- 
tor (FET) structure. When SiC^/Si substrate is used 
as the back-gate insulator of the FET, the tuning range 
of the Fermi level of graphene by the application of the 
gate voltage is at the extent of several ±0.1 eV 18 which 
would be sufficient to alter the chemical reactivity on 
the surface. Besides, the additional charge and the gra- 
dient of electric potential generated by the gate electric 




field are expected to change the polarization of adsorbed 
molecules 19 ' 20 and to modify the charge distribution on 
graphene layers and the adsorbed molecules 21 25 leading 
to, e.g., the change in the ease of migration of molecules 
adsorbed on graphene 26 . Some research 27 29 argues that 
the change in the Fermi level caused by a gate electric 
field activates electrochemical redox reactions and the 
accompanying charge transfer causes hysteresis of the 
source-drain current in graphene FET, yet there has been 
no investigation that elucidates the relation between the 
kinetics of adsorption to graphene and gate electric field. 
In this study we investigated gate-tuned molecular oxy- 
gen adsorption through systematic measurements of con- 
ductivity using mainly bilayer graphene (BLG). 

The back-gated BLG-FETs were fabricated on Si0 2 
(300 nm thick) on heavily n-doped silicon substrate by 
means of photolithography. The channel length and 
width were 6 /im and 3.5 /im, respectively. Prior to mea- 
surement, we repeated vacuum annealing (210°C, 10 h) 
to remove the adsorbed moisture and contaminants on 
the surface until no more changes in the gate-dependent 
conductivity a were eventually seen. After the anneal- 
ing, the BLG-FET exhibited its pristine nature, that is, 
ambipolar transport properties with a conductivity min- 
imum around V g = ^c NP < 8 V giving the " charge neu- 
trality point (CNP)" with electrons and holes in BLG 
being equal in density. The Drude mobility ji{n) was 
estimated to be ~ 1 x 10 3 cm 2 V _1 s _1 from the equa- 
tion n(n) = <r/e|n|, where n is the carrier density with 
n = ( c g/ e )( V g ~ ^cnp) (cg/e = 7 x 10 10 cm^V" 1 ; 
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FIG. 1: (a) Schematic of the measurement cycle. First, the field effect transistor channeled by bilayer graphene (BLG-FET) is 
exposed to gaseous O2 while the gate voltage V^, a d is applied (left panel). Then the system is evacuated and the source-drain 
conductivity of the BLG-FET is measured by sweeping the gate voltage V g (right panel). Subsequently gaseous O2 is again 
introduced, and the cycle is repeated. The gas introduction and evacuation are completed in a shorter time than ~ 10 s to 
prevent additional gas from adsorption. The whole cycle is executed at room temperature, (b) Change of the field effect 
behavior due to the O2 exposure with V gj &d = V (run 1). The gate voltage giving the minimum conductivity (the charge 
neutrality point), is shifted from V g = Vcnp (< 8 V) (marked by a brown triangle, before O2 exposure) to the positive direction 
(green triangles) upon O2 exposure. Circles on the curves represent the conductivity at the hole density of 2.5 x 10 12 cm -2 
that are used to calculate Drude conductivity shown in Figure 4a. (c) The same measurement as in the panel (b) with applying 
the finite V^, a d • All the curves are shifted by —Vcnp (Vcnp = 13, 9, and 11 V for the run of V^, a d = +80, +40, and —50 V, 
respectively) in V g direction, i.e., the charge neutrality points of the pristine graphene without the adsorbed oxygen are taken as 
zero gate voltage. The top panel of (c) represents the a vs V g — Vcnp f° r the pristine graphene. The changes of a vs V g — Vcnp 
curve after a single and a double exposure to O2 (the time duration of a single exposure is 30 s) are shown in the center and 
the lower panel of (c), respectively. Filled circles indicate V s hiit — Vcnp — Vcnp (the shift of the CNP) for each curve. 



c g is the capacitance per unit area for the back-gated 
graphene FET on 300 nm-thick Si02 ) 22 , It is in the 
range of the values for BLG-FET in two-probe config- 
uration previously reported 6 , and therefore we confirm 
that the graphene of the present BLG-FET has few de- 
fects that may extremely enhance the chemical reactivity 
of graphene 30 . 

Next we exposed the BLG-FET to 1 atm of high- 
purity (> 99. 9995%) oxygen in the measurement cham- 
ber at room temperature. Instead of measuring conduc- 
tivity with graphene kept in the O2 environment, we 
performed the short-time interval exposure-evacuation 
cycles schematically shown in Figure la; O2 exposure 
was done under the dc gate voltage V^ a d, followed by 
rapid evacuation in less than 10 s (the physisorbed O2 
molecules would be removed immediately without charge 
transfer), and eventually the a vs V g measurement was 



done with sweeping Vg under vacuum. In this cycle, we 
can rule out the possibility of the additional oxygen ad- 
sorption during sweeping gate voltage for the a vs V g 
measurement since the system was evacuated then. In 
addition, we found that the a vs V g curve did not vary 
under vacuum at room temperature at least for more 
than several hours, so that we can also rule out the pos- 
sibility of the oxygen desorption during a vs V g mea- 
surement (taking ca. 10 min to obtain a single a vs 
Vg curve). Therefore, just repeating the cycles substan- 
tially realizes the long-time O2 exposure under V g ^a, the 
length of which is denoted by total O2 exposure time, 
t. Figure lb shows the change in a vs V g by repeat- 
ing the O2 exposure-evacuation cycles without applying 
gate voltage during O2 exposure (run 1, Vg, a d = V). 
The shift of the charge neutrality point by the amount 
of V s hih(t) = Vcnp(^) — ^cnp toward the positive di- 
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FIG. 2: Time dependence of V^hift and doping density with doping due to the O2 exposure under application of various 
Vg,&d- Solid and dotted curves are the fits based on the H kinetics and the P kinetics (see text), respectively. The curve fitting 
is made in the range of exposure time below 2000 min. The electron-dominant region, where the Fermi level of BLG is higher 
than the CNP (V^ad > Vcnp), is painted blue and the hole-dominant region (Vg, a d < Vcnp) is painted white. Carrier type is 
inverted between electron and hole at the point indicated by arrows (numbers aside correspond to the run number) during the 
evolution of the oxygen adsorption. 



rection was observed (Figure lb), which represents hole 
doping to graphene. Prolonged exposure brought further 
hole doping, and eventually the doping density (induced 
charge by the oxygen adsorption) n^(t) = (c g /e)V s hnt(t) 
reached more than 5 x 10 12 cm -2 within a time scale 
of 10 3 min. Note that any hysteresis as observed in 
graphene FET in moist atmosphere 29 ' 31 34 was not found 
in the observed a vsV g curve, so that we can uniquely de- 
termine Vcnp(^) as a function of t. Another remarkable 
feature is the hole conductivity in highly doped regime, 
V g - Vcnp W < -40 V (i.e., \n\ > 3 x 10 12 cm" 2 ). The 
a vs V g curve distorted and exhibited the sublinear de- 
pendence in this regime for the pristine BLG-FET as can 
be seen in Figure lb. Yet it disappeared only after the 
exposure to O2 for 1 min, whereas the carrier doping has 
not proceeded much at that time. Thus this rapid change 
in conductivity feature is discriminated from the slower 
change causing the shift of the CNP; one possibility is 
that the former is due to rapid oxidation of the metal- 
graphene interface 35 ' 36 , which does not shift the Fermi 
level of graphene but asymmetrically varies the conduc- 
tivity. 

Both Vcnp and the mobility of the BLG-FET with 
the oxygen adsorbed were reset to the value for the pris- 
tine BLG-FET by annealing the 2 -exposed BLG-FET 
in vacuum at 200 °C, indicating that O2 desorption read- 
ily proceeds at high temperature without making any 
defects. By virtue of this reversibility of the oxygen ad- 



sorption, we can repeat the conductivity measurements 
in the exposure-evacuation cycles as described above for 
the same device and compare the results. We addition- 
ally carried out four consecutive measurements under the 
same condition except that a finite gate voltage ^, a d was 
applied during O2 exposure; V^ a d was +80 V (run 2 and 
run 3), —50 V (run 4) and +40 V (run 5). Figure lc 
represents the change of the gate-dependent conductiv- 
ity a vs V g — Vq NP at the initial step of the runs 2, 4, 
5: before O2 exposure, (i.e., after vacuum annealing) 
and after the first and second exposure-evacuation cy- 
cles (the duration time for O2 exposure in each cycle is 
30 s, i.e., t = 0.5, 1 min, after the first and second cycle, 
respectively). All the a vs V g — Vcnp curves collapsed 
onto almost the identical curve in the pristine graphene 
as shown in the top panel of Figure lc. Since Vcnp was 
within 10 + 3 V for each run (see the caption of Figure 1), 
the BLG-FET was realized to exhibit its pristine feature 
before O2 exposure. After the O2 exposure (the center 
and the bottom panel of Figure lc) the gate-dependent 
conductivity changes similarly as was also observed for 
run 1 (Figure lb), but the effect of applying V^ a d during 
O2 exposure is marked by the clear difference in V s hnt{t). 
There is a tendency that Vghift(^) is larger for higher V^ ?a H. 
and smaller for lower V^ a d, indicating that hole doping 
proceeds more intensively to graphene with higher Fermi 
level. This trend is pronounced on the increase in the ex- 
posure time, as confirmed by comparing the center and 
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FIG. 3: Double logarithmic plot of time dependence of 
dn^/dt. Symbols are taken in common with Figure 2. We 
estimate dn^/dt from the differential between the neighbor- 
ing data points for each run in Figure 2. Lines are linear fits 
for all the differential data. The slope of each curve gives u = 
1.02, 0.91, 0.86, and 0.77 (dn^/dt oc t~ u ) for the gate voltage 
V^ad = +80, +40, 0, and -50 V, respectively. 



the bottom panel of Figure 1(c). 

We tracked the temporal evolution of the gate- 
dependent conductivity over a wide time range between 
10°-10 3 min. Figure 2 represents T4hift for runs 1-5 with 
respect to O2 exposure time, in which the correspond- 
ing doping density owing to the oxygen adsorption, n^, 
is also shown in the right axis. The tendency that the 
high Vjx,ad leads to rapid doping can be seen clearly over 
a whole time range; e.g., to reach the doping level of 
K>hift = 40 V, it took ca. 300 min for nonbiased BLG- 
FET. In contrast, hole doping is so enhanced for the 
BLG-FET of V gM = 80 V that it took only 4 min, on the 
other hand so suppressed for that of Vg ?a d = — 50 V that 
it took more than 1000 min. The doping density increases 
almost linearly with respect to logt for V^ a d = +80 V 
and +40 V, whereas superlinearly for lg, a d = V and 
—50 V. The plots for runs 2 and 3, having common 
^g,ad = +80 V, are completely on the same line, which 
verifies that the thermal annealing in vacuum for the re- 
producing of the undoped state in the BLG does not 
affect the behavior of adsorption. 

Figure 3 shows the time dependence of the doping 
rate dn^/dt estimated from the differential between 
the neighboring data points in Figure 2. It is obvi- 
ous that the doping rate changes in accordance with 
dn^/dt oc t~ u . The power u is dependent on V^, a d5 u ~ 1 
for VJr 5 ad = +80 V and it decreases for the runs with 
lower V^ a d. This deviates from the conventional Lang- 
murian kinetics for molecular adsorption which would 
give dn^/dt oc exp(— t/r) with a constant r. 

Careful verification is necessary to inquire the gate- 
voltage-dependent and non-Langmurian temporal change 



of the molecular doping since the rate for doping density 
dn^/dt is related to both of the rate for the chemisorp- 
tion of molecules (d7V ox /d£, where 7V ox is the areal density 
of the adsorbed oxygen molecules) and the transferred 
charge per adsorbed molecule (the charge/molecular ra- 
tio, Z). Therefore, we analyze the mobility that in- 
cludes the information of the scattering mechanism of 
the conducting electrons and the charge of the adsorbed 
molecules. Within a standard Boltzmann approach 37 , 
the mobility is changed inversely proportional to the 
density of the scattering centers (i.e., the adsorbed 
molecules), 7V ox . In the realistic case, the inverse mobil- 
ity is given as a function of N ox and the carrier density 
n, which reads 5 

fJL(n,N m ) C(n) /i (n) 

Here /io(^) represents the mobility of the pristine 
graphene without the adsorbed oxygen. The coefficient 
C(n) represents the feature of carrier scattering by the 
adsorbed oxygen. On the one hand, the charged-impurity 
scattering 38 gives C(n) oc [1 + 6.53^ (d + X T f)]/Z 2 
for the BLG in the low-carrier-density regime and 
in the limit of d -+ within the Thomas-Fermi 
approximation 38 , where d is the distance between the 
impurities and the center of the two layers of the BLG 
(see the inset of Figure 4a for the definition), and the 
screening length Atf = fth 2 /4m*e w 1 nm (k: dielectric 
constant) 38 . On the other hand, the short-range delta- 
correlated scatterers give the constant C(n) = C s 38 , or 
the strong impurities with the potential radius R give 
C(n) oc \\n(RyJim)] 2 in the off-resonance condition 39 , 
which is a decreasing function of n in the regime of 
n ~ 10 12 cm -2 taking R to be several angstroms. As for 
the relation between the concentration of the adsorbed 
oxygen molecules 7V ox and the 02-induced doping den- 
sity n^x", we assume that the charge Ze of each adsorbed 
oxygen molecule dopes the carriers — Ze in the BLG, no- 
tably, tt^x" = —ZN x. To be exact, the amount of the 
induced charge is not such a simple function 6 propor- 
tional to the number of adsorbed molecules due to the 
energy-dependent DOS of BLG 40,41 and the anomalis- 
tic screening effect therein 5 ' 6,42 . Yet in the low energy 
regime of BLG where the DOS is envisaged to be con- 
stant, the assumption above is appropriate. 

Figure 4a shows the inverse Drude mobility vs 
n^x" plots at the carrier density of n = 2.5 x 10 12 cm -2 
(marked by the open circles in Figure lb for run l:V gya a = 
V) for V^ad = +80, +40, 0, and -50 V. Linear increase 
in with respect to was found. This, along with 
the linearity between and 7V ox given by Eq. (1), im- 
plies that Z is not a function of N ox , i.e., invariant against 
the increase of the adsorbed molecules. Interestingly, the 
slope of vs tt^x" plot (the inverse of the slope corre- 
sponds to C(n)\Z\ in Eq. (1)) depends on V^ a d- Note 
that before O2 was introduced (n^ = 0), we observed 
& 28 Vsm -2 irrespective of V^ a d, and thus the dif- 
ference in the mobility by V^ >a d genuinely results from 
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FIG. 4: (a) Inverse mobility vs the doping density 
for each run in Figure 2. Symbols are taken in common with 
Figure 2. Lines are linear fits (as for V gj &d = +80 V, the 
data for both run 2 and run 3 are included). Inset: Adsorbed 
oxygen species with charge Ze positioned at the distance d 
away from the center of BLG. (b) C\Z\ vs the carrier density 
n (see Eq. (1)). For V g , ad = +80 V, C\Z\ is acquired by 
gathering the data for run 2 and run 3. Theoretical results 
based on a charged- impurity scattering mechanism are shown 
with the various distance d and charge/molecular ratio Z: the 
solid line represents the result for d — 0.43 nm and Z — 0.38, 
the dashed line for d = 0.73 nm and Z — 0.38, and the dotted 
line for d — 0.43 nm and Z — 0.28, respectively. 



the adsorbed oxygen instead of other unintentional im- 
purities on the BLG or the Si02 substrate. In Figure 4b, 
the inverse of the slope, C(n)\Z\ 1 is shown for the var- 
ious carrier densities, n. Therein we omit the data in 
the low carrier regime of n < 2.5 x 10 12 cm -2 , in which 
the residual carriers due to electron-hole puddles can- 
not be disregarded and the carrier density n (and thus 
also the Drude mobility) cannot be correctly estimated 
only by considering the gate electric field effect 43 . Sim- 
ilarly to the charged impurity model rather than oth- 
erwise, C(n)\Z\ is increasing with n. The dependence 
experimentally observed, however, still deviates from the 
theoretical calculated results within the charged impurity 
model plotted in Figure 4b for various d and Z (assuming 
d and Z are invariant to n). 



The difference in C(n)\Z\ depending upon Vg ja d indi- 
cates that the electronic polarity of graphene varies the 
adsorption states of oxygen molecules, leading to the vari- 
ation of d and Z. When the positive (negative) V^, a d is 
applied, negative (positive) carrier is electrically induced 
on graphene, which may modify the interaction between 
graphene and the adsorbed oxygen molecules with the 
negative charge, e.g., the Coulomb interaction and the 
overlap of the orbit als. Eventually, the stable adsorption 
state is varied by V^ >a d, leading to the difference in mo- 
bility. Besides, let us recall that conductivity measure- 
ment process is set apart from the O2 adsorption process 
and that constant ^, a d is not applied when the mobility 
is measured (Figure la). Accordingly, whereas the sta- 
ble adsorption state of oxygen during the conductivity 
measurement may differ from that during adsorption, the 
adsorbed oxygen molecules are kept in the former state 
during conductivity measurement, and the mobility vary- 
ing by Vg s ad is actually observed. This indicates that the 
energetic barrier exists for charge redistribution between 
graphene and the adsorbed oxygen molecules (shown be- 
low), and once the adsorption is accomplished, the charge 
Ze on each adsorbed oxygen molecules will not immedi- 
ately change just after switching on/off the gate voltage. 
One possible reason for the deviation between the ex- 
perimental results and theoretical curve is that d varies 
accompanied with the change in n (or sweeping V g ). Yet 
actually, since modifying d by several angstroms results 
in the small change in (C(n)\Z\) as shown in Figure 4b, 
it is necessary to investigate more about the behavior of 
the adsorbed oxygen molecules in the gate electric field 
in a the future study. 

It is controversial what kind of oxygen species does ac- 
tually cause the hole doping to graphene. Because the 
electron affinity of O2 (0.44 eV 44 ) is much lower than 
the work function of graphene (4.6 eV 45 ), direct charge 
transfer between them seems unfavorable. Instead, with 
an analogy of the charge doping of diamond surface 46 , 
there is a widely accepted 28,33,47 hypothesis that the hole 
doping proceeds through an electrochemical reaction 46 
such as: O2 + 2H20 + 4e~ = 40H~, by which the charge 
transfer is favorable due to the lowered free energy change 
AG = —0.7 eV 28 on the condition that the oxygen pres- 
sure is 1 atm and pH = 7. This electrochemical reaction 
needs the aid of water that is mostly eliminated in the ex- 
periment by annealing (we observed no hysteresis in the 
g vs V g curve, that is, there are few charge traps often 
attributed to residual moisture on graphene or its sub- 
strate). Yet as for graphene deposited on the hydrophilic 
Si02 substrate 32 , it is possible that a small amount of 
residual water molecules (more than the chemical equiv- 
alent of O2) are trapped on Si02 surface or voids, which 
cannot be easily removed by the vacuum annealing at 
200 °C in comparison with those on graphene surface. 
We suggest that the electrochemical mechanism is plau- 
sible also in our case, yet the adsorbed molecule could 
be other chemical species than OH - , the charge of which 
may be dependent on V^ a d- 
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FIG. 5: Schematic energy diagrams of the kinetics of O2 
adsorption (H kinetics), (a) Path for electron transfer in 
this model is shown by the blue dotted arrow; electrons in 
BLG (the electrochemical potential Cg) are transferred to 
O2 molecules via the transition state (the circled T at the 
level of Cts ), giving the adsorbed oxygen species (the cir- 
cled A at the level of Cads)- The activation energy, the free 
energy change, and the level of the CNP are denoted by *E, 
AG and Ccnp, respectively. The Fermi level is defined by 
£f = (g — Ccnp- As a demonstration, the case for Vg, a d > is 
presented, (b) Temporal change in the activation energy and 
the Fermi energy due to the adsorption of oxygen molecules 
to BLG negatively doped by the positive gate voltage (as in 
case of panel a). The left panel represents the case before 
oxygen adsorption (t = 0) , and the right panel represents the 
oxygen exposure for the time t. Here the energy is measured 
from the CNP. The red and the blue tick marks denote the 
level for the transition state and the Fermi level, respectively. 
Oxygen adsorption lowers the Fermi level, accompanied with 
the increase in hole doping of (equal to the area of grayed 
part, fs*$ D(sf) d£F), and the activation energy increases 
according to Eq. (2) 



In light of the discussion above, the electrochemical 
description 33 ' 48 is expected to be applicable for the ob- 
served adsorption kinetics of oxygen to the BLG. Here 
we premise that the molecular adsorption is determined 
by the electrochemical potential of graphene, and con- 
sider the charge transfer kinetics in an approach based 
on Butler- Volmer theory 49,50 . The model is schemati- 
cally depicted in Figure 5. The probability of the adsorp- 
tion reaction is determined by the electrochemical poten- 
tial of graphene (Cg) and that in the equilibrium condi- 
tion of the oxygen-chemisorption reaction (Cads)- When 



AG = Cads — Cg < 0, the electrons favorably transfer 
from of graphene to the adsorbed oxygen (denoted by 
"A" in Figure 5a), and the oxygen-adsorption reaction 
proceeds. For charge transfer, the electrons should go 
through some energy barrier; we assume that electrons 
tunnel from BLG to the O2 molecules via a single transi- 
tion state (denoted by "T"), whose electrochemical po- 
tential is Cts- The difference $E = Cts — Cg corresponds 
to the activation energy of the oxygen-chemisorption re- 
action, which determines the frequency of the electron 
transfer. Whereas Cg is dependent on the Fermi level e F 
as Cg = £f + Ccnp (Ccnp is the electrochemical potential 
of the CNP), we envisage that Cts (or is a function 
of £y as well. In the framework of the Butler- Volmer 
theory, we obtain the dependence of $E on Ev as 



E(e F + de F ) = *E(ev) - ol de^ 



(2) 



where a(> 0) is a constant related to the "transfer coef- 
ficient" in the Butler- Volmer theory that associates the 
activation energy with the electrochemical potential (not 
the Fermi energy) ; thus herein we call a as "pseudo trans- 
fer coefficient" (see Supporting Information for detail in 
the derivation). That is, we have the assumption that 
the activation energy scales linearly with the Fermi en- 
ergy. Further assuming that the molecular adsorption 
rate d7V ox /dt is controlled by the electron transfer pro- 
cess and is not strongly affected by other contributions 
such as molecular diffusion 14 , it is given by 



dN , 
dt 



= X D (e P + t£(e P )) / (e P + *£(£ F ); e P ) 



(3) 



where D(e) is the density of states (DOS) of BLG (e is 
the energy measured from Ccnp) and f(e;ev) = [1 + 
(e — £y) / k^T)~ l is the Fermi-Dirac distribution func- 
tion. The coefficient x does not depend on e F (if the 
distance between the adsorbing molecules and BLG var- 
ied depending on e F or lg, a d 7 the tunneling frequency 
would be affected so that \ might be dependent on them 
as well; yet herein we ignore such effect for simplicity). 
The right-side of Eq. (3) represents the tunneling rate of 
the electron from the graphene to oxygen at the energy 
level of transition state, e = e F + ^E(e F ) = Cts — Ccnp- 
According to Eq. (2) and Eq. (3), the molecular adsorp- 
tion rate is dependent on the Fermi level of graphene. 
Using them, we can explain both the temporal evolution 
of the doping rate and its dependence on the gate voltage 

^g,ad- 

Let us discuss the temporal change in the doping rate. 
The Fermi level of graphene is lowered with the increase 
of the adsorbed oxygen molecules, because the positive 
charge is induced on graphene by the charge Ze they 
possess. Recalling that 77^ = — Z7V ox , the doping rate 
is given by dn^/dt. Furthermore, since $E ^> k^T is 
fulfilled as shown later, we also approximate that f(ev + 
$E(eF);s F ) ~ ex.p(—$E(s F )/kBT). Using the relation 
dn^ = — D(£F)d£F, we acquire a formula describing the 
temporal change of the Fermi level: 
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- D Mt)) = (e F (t) + *E (e F (t))) exp ( ^~^ ) , (4) 



where £f(£) = £f(£, ^,ad), expressing that the Fermi 
level is a function of the exposure time t and the gate 
voltage ^g,ad 7 and £f(0) = £F(0,Vg, a d) 7 the Fermi level 
at t — 0. We have specifically defined two constants, the 
pseudo transfer coefficient a te (the subscript "te" abbre- 
viates "temporal evolution" ) and 

P = -^F^\-^T-) (5) 

The right side of Eq. (4) represents the product of the 
charge transfer frequency and the amount of charge per 
adsorbed molecule, whereas the left side does the resul- 
tant amount of the doped charge. Because BLG (or also 
single layer graphene) has the low DOS around the CNP 
compared to metal, the small amount of carrier doping 
results in the large shift in the Fermi level, which effec- 
tively controls the kinetics. Thus the adsorption kinetics 
is well described by Eq. (4), the equation focusing on the 
Fermi level. When we envisage that the BLG is approx- 
imately described by two-dimensional parabolic disper- 
sion of the free electron, the DOS becomes constant as 
D = Dp = 7±/7t(^f) 2 , where vp = (V3/2)joa/h is the 
Fermi velocity in SLG, a = 2.46 A is the in-plane lattice 
constant, and 70 = 3.16 eV and 7^ « 0.4 eV 51,52 are the 
intrasheet and inter sheet transfer integrals, respectively. 
In this case (hereafter labeled as P kinetics), Eq. (6) is 
readily integrated, giving 

n^{t) = ~{s F (t) - £f(0)}D p = l n (i +pt ) 

\ "te / 

(6) 

Note that Eq. (6) satisfies dn ox /dt oc 
exp[— at e n^i(t) / k^T D) and is equivalent to the in- 
tegrated form of Elovich equation 53,54 , the empirical 
equation that is widely applicable to chemisorptions 
onto semiconductors. When the hyperbolic DOS of 
BLG 41 is reflected to Eq. (4), more accurate but more 
complicated expression of Ep(t) is acquired (denoted as 
H kinetics), given by 

" S eXP ( ate ^T ) S ^' £F ^ +P t = > (7) 

where S[£p(0), sp(t)] is a function that depends on the 
DOS at the Fermi level and that at the transition state 
level (derived in the Supporting Information). 

We performed curve fitting of the experimental results 
of n ox (t) with Eq. (6) and Eq. (7) for P kinetics and 
H kinetics, respectively. The difference in V^ a d by runs 
is simulated by the dependence of £f(0) on V^ a d first 
without considering the gap-opening effect 55 ' 56 , that is, 



I 

we calculate sp(0) using the relation that the charge of 
Cg(V^ a d — ^c NP ) doped on BLG by applying V^ a d is equal 

to Jq F( ^ D(e)de. Irrespective of the kinetic models, All 
the theoretical curves are well fitted to the experimental 
results except those for run 4 in the range above 10 3 
min. Eq. (4) is invalid in the first place in the long- 
time regime in which the adsorption rate is almost as 
low as the desorption rate, since the present treatment 
includes no contribution of desorption. The deviation 
is, however, contrary to the expectation; the desorption 
should suppress the evolution of the hole doping yet the 
enhanced doping was actually observed. Thus we suspect 
that it is due to long-time scale chemisorption of oxygen 
onto graphene which is ubiquitously observed in carbon 
materials 57 ' 58 . Note that the several volts of offsets are 
added in T^hift (corresponding to the doping density of 2 x 
10 11 cm -2 ) to improve the fitting in the time range of t < 
10° min. This small offset corresponds to a fast reaction 
that finishes at the very initial stage of adsorption, e.g., 
due to the reactive chemisorption of O2 to the defect or 
edge site of graphene 59 ' 60 . 

Figure 6 shows the profiles of the fitting results. Prac- 
tically, the fitting parameters are following two: a te and 
p, but in order to derive the initial activation energy be- 
fore the oxygen adsorption, $E(ep(0)), from p based on 
Eq. (6), we briefly assume that the charge/molecular ra- 
tio is independent of the gate voltage and determine that 
Zxi± = 1 x 10 5 eV 2 min -1 from the pre-exponential fac- 
tor in the literature 33 , c ox ^ei ~ 10 17 cm _2 s _1 (c 
in the literature corresponds to X7±/ 7F (^f) 2 in this pa- 
per). Note that if we assume another value than Zx^± = 
1 x 10 5 eV 2 min _1 , it results in a uniform shift of the 
calculated $E(ep(0)). In addition, the charge/molecular 
ratio of the adsorbed molecule Z is likely dependent on 
V gydL d from the discussion about the mobility, but it leads 
the shift of tE(e F (0)) only by - k B T = 0.026 eV. On 
the one hand, a te (representing the temporal change of 
the activation energy) exhibits a significant deviation be- 
tween the H kinetics and the P kinetics (Figure 6a), or de- 
pends on the treatment of the DOS of BLG. This devia- 
tion is understood as follows: H kinetics reflects the DOS 
of BLG that is a monotonically increasing function with 
respect to \ep\ having a minimum (D = Dp) at the CNP 
(Figure 6(c)), yet P kinetics does not. Since the down- 
ward shift of the Fermi level upon the oxygen adsorption 
depends on the DOS at the Fermi level, the range of the 
change in the Fermi level differs between two kinetics 
models (Figure 6(c)). Thus a te , related with the Fermi 
level by Eq. (2), is calculated differently. This result is 
contrasted with the conventional electrochemical reaction 
on the metal electrodes in which the kinetics is not signif- 
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FIG. 6: Summary of parameters obtained by the curve fitting 
in Figure 2 with P kinetics (P) and H kinetics (H) compared, 
(a) Dependence of the pseudo transfer coefficient ate on lg, a d- 
The pseudo transfer coefficient a g f acquired by the linear fit- 
ting to panel b is also plotted, (b) Initial activation energy 
*E(ef(0)) plotted as a function of ef(0) — £f(0, V g ). (c) (Left 
panel) Density of states (DOS, solid line) for BLG as a func- 
tion of the Fermi energy, where Dp (DOS at the energy of 
£f = 0) is taken as unit of the DOS. The solid line denotes 
the DOS applied for the H kinetics. In P kinetics the DOS 
is treated as constant; D = Dp (dotted line). (Right panel) 
Change of the Fermi level by oxygen adsorption as a function 
of V^ad- Bars and boxes represent the range calculated under 
P kinetics and H kinetics, respectively, assuming 7^ = 0.4 eV. 
The Fermi level before the oxygen adsorption (sf(0)) is at the 
point indicated by triangles (black squares) for H (P) kinetics, 
whereas the Fermi level shifts downward with the increase of 
the adsorbed oxygen, and eventually reaches the bottom end 
of the boxes (bars) in the end of the oxygen adsorption for 
the H (P) kinetics. The energy is measured from the CNP 
(dashed line), and the gap-opening effect is not considered. 



icantly affected by the DOS of the electrodes and is just 
owing to the low DOS of the BLG with comparison to the 
metal. On the other hand, we found that ^E(e^(0)) is a 
decreasing function of £f(0) = £f(0, V^ a d) (shown in Fig- 
ure 6b, therein we plot ^E(sf(0)) with respect to £f(0) 
instead of V^ ?a d)- From Eq. (2), the slope of the plots in 
Figure 6b corresponds to the pseudo transfer coefficient 
a, and we found that a g f = 0.36 and 0.42 for P kinetics 
and H kinetics, respectively (the subscript "gf" abbrevi- 
ates "gate electric field" ) . Herein we distinguish a g f from 
ate', being different from a te that is calculated based on 



the temporal Fermi level shift by oxygen adsorption (thus 
ate inevitably includes the oxygen adsorption effect), a g f 
is acquired by tuning the Fermi level electrically at t = 
before oxygen adsorption. We found that a g f is much 
smaller than at e (> 1), or rather near 0.5, a typical trans- 
fer coefficient for simple redox reactions 48 ' 49 (Figure 6a). 
Specifically, the oxygen adsorption effect included only 
in ate but not in a g f is attributed, e.g., to the electric 
dipole layer 61 formed between graphene and the adsorbed 
molecules and the Coulomb interaction between the ad- 
sorbed oxygen molecules. We expect that these effects 
also raise the activation energy roughly in proportional 
to the number of molecules N ox (or the doping density 
n^), and thus we have the expression for the additional 
adsorption effect: d*E = £ mo idn^ = — £ m o\D(eF)deF 
(£moi* the proportional coefficient; herein we use the 
relation dn^ = —D(sF)dsv again). Then we acquire 
at e — Cigf + €moi(D(eF)) ((D(ev)) denotes the average 
DOS in the range of the Fermi level for each run), which 
indicates that a te is large when the DOS is large in 
the level far away from the CNP (see Figure 6c). In- 
deed, as shown in Figure 6a, a te for H kinetics shows 
the V-shaped dependence on Vjr ?a d with the minimum for 
^g,ad = +40 V, which gives the smallest average DOS. 
Though the kinetics for the molecular adsorption is af- 
fected by various effects as mentioned above, we represent 
them by the parameter a te and succeed in accounting for 
the observed kinetics in a facile way. 

Finally let us account for the power-law dependence 
of dn^/dt ex t~ u shown in Figure 3. It is helpful to 
look on the simpler P kinetics for the assessment of 
dn^/dt; from Eq. (6), we obtain dn^/dt ex p/(l +pt). 
Approximately we have dn^/dt ex (ii(t) n-1 )t -n where 
u ~ [1 + < 1 (the time (t) is the center of 

the expansion of ln(dn^/d£) in terms of lnt, and it is a 
good approximation if p(t) ^> 1, or otherwise if p(t) ~ 1 
in the time range such that t = [10 — 1 (t), 10(t)]). Since 
p is intensively dependent on V^ a d (recall that (i) p is 
exponentially decaying against *E(ev(0)) as represented 
in Eq. (5), (ii) *E(ev(0)) linearly decreases with respect 
to £f(0) with the slope of — a g f as shown in Figure 6c, 
(hi) £f(0) is an increasing function of V^ a d. We acquired 
p =19.6, 2.9, 0.98 and 0.32 for V gM = +80, +40, 0, 
and —50 V, respectively, by fitting within the P kinetics 
model), we can find that u is almost unity for a positively 
high V^ad and tends be smaller for negatively high V^ ?a d 
within the time range experimentally scoped, which is 
consistent with the observed behavior. When the elec- 
trochemical mechanism governs the kinetics of the oxy- 
gen adsorption, the activation energy of the charge trans- 
fer continually increases with the O2 exposure time in- 
creasing. This effect leads to non-Langmurian kinetics 
of the oxygen adsorption and the power-law decrease of 
dn^/dt, even though neither the desorption process nor 
the saturation limit of adsorption are taken into consid- 
eration. 

In BLG, it is known that a band gap opens due to the 
energy difference between two layers 9 . The gate electric 
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field as well as the adsorbed oxygen may produce such 
a strong energy difference that the eventual band gap 
should affect the time evolution of molecular adsorption. 
The band gap opening effect is expected to exhibit most 
prominently when the Fermi level goes across the CNP 
(at the point shown by arrows in Figure 2), while we can- 
not find such behavior obviously. We guess it is partly be- 
cause most of the adsorbed oxygen molecules exist in the 
interface between graphene and the gate dielectric Si02- 
For the band gap opening effect to appear, it is neces- 
sary that the gate electric field and the molecular field 
enhance each other when the Fermi level is near the CNP 
(i.e., the charge induced by the gate electric field and that 
by the adsorbed molecules including unintentional resid- 
ual impurities on the Si02 substrate are balanced), yet it 
is possible only when the molecules mainly adsorb on the 
top surface of BLG, and not for the molecules adsorbed 
in the interface (Figure S2, Supporting Information). Or 
it may be partly because the disordered potential due 
to the impurities in the substrate fluctuates the energy 
level around which the band gap exists 62 , eventually blur- 
ring band gap opening effects and chemical reactivity 63 
of BLG. Details about the band-gap opening effects are 
discussed in the Supporting Information. 

In summary, we investigated the weak chemisorption of 
O2 molecules on bilayer graphene by measuring its trans- 
port properties. The hole doping due to O2 chemisorp- 
tion is remarkably dependent on the gate voltage, and 
the amount of the doped carrier increases with O2 ex- 
posure time, the rate of which is in accordance with 
oc t~ u (u < 1) rather than with conventional Langmuirian 
kinetics. We conclude from these that an electrochem- 



ical reaction governs the O2 chemisorption process, in 
which the rate of the chemisorption is determined by 
the Fermi level of graphene, and indeed succeed in ac- 
counting for the observed kinetics by the analysis based 
on the Butler- Volmer theory. We also found that the 
chemisorbed molecules decrease the mobility of graphene, 
and interestingly, the mobility change is dependent on 
the gate voltage applied during the adsorption, indicat- 
ing that the adsorption state, e.g., transferred charge or 
distance between a molecule and graphene, can be mod- 
ified electrically. Graphene, offering a continuously tun- 
able platform for study of chemisorptions on it, realizes 
the electrical control of the adsorption by gate electric 
field, a novel and versatile method in which we would 
explore extensively a wider variety of host-guest interac- 
tions between graphene and foreign molecules. 

Supporting Information. Additional descrip- 
tions about (i) the experimental method, (ii) the 
electrochemistry-based kinetics model and the mathe- 
matical derivation for it, and (iii) an expanded discussion 
about the gap-opening effects. 

Acknowledgments 

The authors acknowledge support from Grant-in- Aid 
for Scientific Research No. 20001006 from the Ministry 
of Education, Culture, Sports, Science and Technology, 
Japan. The authors thank M. Kiguchi, T. Kawakami, K. 
Yokota, and Y. Kudo for useful discussion. 



* Electronic address: mail to: sato.y.an@rn.titech.ac.jp 

1 R. Balog, B. J0rgensen, L. Nilsson, M. Andersen, 
E. Rienks, M. Bianchi, M. Fanetti, E. Laegsgaard, 

A. Baraldi, S. Lizzit, Z. Sljivancanin, F. Besenbacher, 

B. Hammer, T. G. Pedersen, P. Hofmann, and 
L. Hornekaer, Nat. Mater. 9, 315 (2010). 

2 Y. Zhang, A. Kolmakov, S. Chretien, H. Metiu, and 
M. Moskovits, Nano Lett. 4, 403 (2004). 

3 A. H. Castro Neto, F. Guinea, N. M. R. 
Peres, K. S. Novoselov, and A. K. Geim, 
Rev. Mod. Phys. 81, 109 (2009). 

4 F. Schedin, A. K. Geim, S. V. Morozov, E. W. Hill, 
P. Blake, M. I. Katsnelson, and K. S. Novoselov, 
Nat. Mater. 6, 652 (2007). 

5 J. H. Chen, C. Jang, S. Adam, M. S. Fuhrer, E. D. 
Williams, and M. Ishigami, Nat. Phys. 4, 377 (2008). 

6 S. Xiao, J.-H. Chen, S. Adam, E. D. Williams, and M. S. 
Fuhrer, Phys. Rev. B 82, 041406 (2010). 

7 J. H. Chen, C. Jang, M. Ishigami, S. Xiao, 
W. G. Cullen, E. D. Williams, and M. S. Fuhrer, 
Sol. Stat. Commun. 149, 1080 (2009). 

8 D. W. Boukhvalov and M. I. Katsnelson, 
Appl. Phys. Lett. 95, 023109 (2009). 

9 E. V. Castro, K. S. Novoselov, S. V. Morozov, 



N. M. R. Peres, J. M. B. Lopes dos Santos, J. Nils- 
son, F. Guinea, A. K. Geim, and A. H. Castro Neto, 
J. Phys.: Condens. Matter 22, 175503 (2010). 

10 D. A. Abanin, A. V. Shytov, and L. S. Levitov, 
Phys. Rev. Lett. 105, 086802 (2010). 

11 R. M. Ribeiro, N. M. R. Peres, J. Coutinho, and P. R. 
Briddon, Phys. Rev. B 78, 075442 (2008). 

12 T. Lohmann, K. von Klitzing, and J. H. Smet, 
Nano Lett. 9, 1973 (2009). 

13 V. V. Cheianov, O. Syljuasen, B. L. Altshuler, and V. I. 
Fal'ko, EPL 89, 56003 (2010). 

14 H. E. Romero, P. Joshi, A. K. Gupta, H. R. Gutier- 
rez, M. W. Cole, S. A. Tadigadapa, and P. C. Eklund, 
Nanotechnology 20, 245501 (2009). 

15 T. O. Wehling, M. I. Katsnelson, and A. I. Lichtenstein, 
Chem. Phys. Lett. 476, 125 (2009). 

16 K. Takai, S. Eto, M. Inaguma, T. Enoki, 
H. Ogata, M. Tokita, and J. Watanabe, 
Phys. Rev. Lett. 98, 017203 (2007). 

17 K. M. McCreary, K. Pi, A. G. Swartz, W. Han, W. Bao, 
C. N. Lau, F. Guinea, M. I. Katsnelson, and R. K. 
Kawakami, Phys. Rev. B 81, 115453 (2010). 

18 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 



10 



A. A. Firsov, Nature 438, 197 (2005). 

19 J. Zhou, Q. Wang, Q. Sun, P. Jena, and X. S. Chen, 
Proc. Natl. Acad. Sci. USA 107, 2801 (2010). 

20 P. Saalfrank, J. Chem. Phys. 113, 3780 (2000). 

21 W. Liu, Y. H. Zhao, J. Nguyen, Y. Li, Q. Jiang, and E. J. 
Lavernia, Carbon 47, 3452 (2009). 

22 Y.-H. Lu, L. Shi, C. Zhang, and Y.-P. Feng, 
Phys. Rev. B 80, 233410 (2009). 

23 V. W. Brar, R. Decker, H.-M. Solowan, Y. Wang, 
L. Maserati, K. T. Chan, H. Lee, c. O. Girit, A. Zettl, 
S. G. Louie, M. L. Cohen, and M. F. Crommie, 
Nat. Phys. 7, 43 (2011). 

24 J.-H. Parq, J. Yu, Y.-K. Kwon, and G. Kim, 
Phys. Rev. B 82, 193406 (2010). 

25 X. Tian, J. Xu, and X. Wang, 
J. Phys. Chem. B 114, 11377 (2010). 

26 A. M. Suarez, L. R. Radovic, E. Bar-Ziv, and J. O. Sofo, 
Phys. Rev. Lett. 106, 146802 (2011). 

27 A. A. Kaverzin, S. M. Strawbridge, A. S. Price, 
F. Withers, A. K. Savchenko, and D. W. Horsell, 
Carbon 49, 3829 (2011). 

28 H. Pinto, R. Jones, J. P. Goss, and P. R. Briddon, 
Phys. Status. Solidi. A 207, 2131 (2010). 

29 Y. Yang and R. Murali, 
Appl. Phys. Lett. 98, 093116 (2011). 

30 P. L. de Andres and J. A. Verges, 
Appl. Phys. Lett. 93, 171915 (2008). 

31 Y. Shi, W. Fang, K. Zhang, W. Zhang, and L.-J. Li, 
Small 5, 2005 (2009). 

32 C. M. Aguirre, P. L. Levesque, M. Paillet, F. La- 
pointe, B. C. St-Antoine, P. Desjardins, and R. Martel, 
Adv. Mater. 21, 3087 (2009). 

33 P. L. Levesque, S. S. Sabri, C. M. Aguirre, J. Guillemette, 
M. Siaj, P. Desjardins, T. Szkopek, and R. Martel, 
Nano Lett. 11, 132 (2010). 

34 S. S. Sabri, P. L. Levesque, C. M. Aguirre, 
J. Guillemette, R. Martel, and T. Szkopek, 
Appl. Phys. Lett. 95, 242104 (2009). 

35 R. Nouchi and K. Tanigaki, 
Appl. Phys. Lett. 96, 253503 (2010). 

36 T. Yamada, Phys. Rev. B 69, 125408 (2004). 

37 S. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, 
Rev. Mod. Phys. 83, 407 (2011). 

38 S. Adam and S. Das Sarma, 
Phys. Rev. B 77, 115436 (2008). 

39 M. Monteverde, C. Ojeda-Aristizabal, R. Weil, 
K. Bennaceur, M. Ferrier, S. Gueron, C. Glat- 
tli, H. Bouchiat, J. N. Fuchs, and D. L. Maslov, 
Phys. Rev. Lett. 104, 126801 (2010), Supplementary 
Material II. 

40 E. McCann and V. I. Fal'ko, 
Phys. Rev. Lett. 96, 086805 (2006). 

41 S. Adam and M. D. Stiles, 



Phys. Rev. B 82, 075423 (2010). 

42 S. Adam, E. H. Hwang, V. M. Galitski, and S. Das Sarma, 
Proc. Natl. Acad. Sci. USA 104, 18392 (2007). 

43 Y. W. Tan, Y. Zhang, K. Bolotin, Y. Zhao, S. Adam, 
E. H. Hwang, S. Das Sarma, H. L. Stormer, and P. Kim, 
Phys. Rev. Lett. 99, 246803 (2007). 

44 K. M. Ervin, I. Anusiewicz, P. Skurski, J. Simons, and 
W. C. Lineberger, J. Phys. Chem. A 107, 8521 (2003). 

45 S. J. Sque, R. Jones, and P. R. Briddon, 
Phys. Status Solidi A 204, 3078 (2007). 

46 V. Chakrapani, J. C. Angus, A. B. Anderson, S. D. 
Wolter, B. R. Stoner, and G. U. Sumanasekera, 
Science 318, 1424 (2007). 

47 S. Ryu, L. Liu, S. Berciaud, Y.-J. Yu, H. Liu, P. Kim, 
G. W. Flynn, and L. E. Brus, Nano Lett. 10, 4944 (2010). 

48 I. Heller, J. Kong, K. A. Williams, C. Dekker, and S. G. 
Lemay, J. Am. Chem. Soc. 128, 7353 (2006). 

49 J. O. M. Bockris and Z. Nagy, 
J. Chem. Educ. 50, 839 (1973). 

50 J. O. Bockris, A. K. N. Reddy, and M. Gamboa-Aldeco, 
Fundamentals of Electrodics, 2nd ed., Modern Electro- 
chemistry, Vol. 2A (Plenum Press, 2000). 

51 W. W. Toy, M. S. Dresselhaus, and G. Dresselhaus, 
Phys. Rev. B 15, 4077 (1977). 



M. 



S. Dresselhaus and G. Dresselhaus, 



Adv. Phys. 51, 1 (2002). 

53 C. Aharoni and F. C. Tompkins, Kinetics of adsorption 
and desorption and the Elovich equation, edited by D. D. 
Eley, H. Pines, and P. B. Weisz, Vol. 21 (Academic Press, 
New York and London, 1970) pp. 1-49. 

54 I. S. Mclintock, Nature 216, 1204 (1967). 

55 T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Roten- 
berg, Science 313, 951 (2006). 

56 E. V. Castro, K. S. Novoselov, S. V. Morozov, 
N. M. R. Peres, J. M. B. L. dos Santos, J. Nilsson, 
F. Guinea, A. K. Geim, and A. H. Castro Neto, 
Phys. Rev. Lett. 99, 216802 (2007). 

57 N. Kobayashi, T. Enoki, C. Ishii, K. Kaneko, and 
M. Endo, J. Chem. Phys. 109, 1983 (1998). 

58 G. U. Sumanasekera, G. Chen, K. Takai, J. Joly, 
N. Kobayashi, T. Enoki, and P. C. Eklund, 
J. Phys. Condens. Matter 22, 334208 (2010). 

59 X. Wang and H. Dai, Nat. Chem. 2, 661 (2010). 

60 F. M. Koehler, A. Jacobsen, K. Ensslin, C. Stampfer, and 
W. J. Stark, Small 6, 1125 (2010). 

61 G. Giovannetti, P. A. Khomyakov, G. Brocks, V. M. 
Karpan, J. van den Brink, and P. J. Kelly, 
Phys. Rev. Lett. 101, 026803 (2008). 

62 G. M. Rutter, S. Jung, N. N. Klimov, D. B. Newell, N. B. 
Zhitenev, and J. A. Stroscio, Nat. Phys. 7, 649 (2011). 

63 R. Sharma, J. H. Baik, C. J. Perera, and M. S. Strano, 
Nano Lett. 10, 398 (2010). 



